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Abstract. The dichotomy between a universal mass function (IMF) and a variable IMF which depends on local physical 
parameters characterises observational and theoretical stellar astronomy. In this contribution the available distributions of 
probability are briefly reviewed. The physical nature of two of them, gamma variate and lognormal, is then explained once 
the framework of the fragmentation is introduced. Interpolating techniques are then applied to the sample of the first 10 pc 
and to the open cluster NGC6649: in both cases lognormal distribution produces the best fit. The three power law function has 
also been investigated and visual comparison with an artificially generated sample of 100000 stars suggests that the variations 
in the spectral index are simply due to the small number of stars available in the observational sample. In order to derive the 
sample of masses, a new formula that allows us to express the mass as a function of the absolute magnitude and (B-V) for 
MAIN V , GIANTS III and SUPERGIANTS I is derived. 
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1. Introduction 

The first analytical expression adopted to explain the low 
mass distribution of stars, see [Salpeter (1955)| was a power 
law of the type ^{M) oc where (,{A4) represents the 

probability of having a mass between and A4 + dA4 ; 
this is called initial mass function, in the following IMF. In 
the range IOMq > M > IMq, the value of a, 2.35, 
has changed little over the decades and a recent evaluation 
quotes 2.3, seeQCroupa (20 01)] A fractional exponential law 
aimed at explaining the mass distribution of asteroids, stars 
and galaxies was then introduced by Kiang (1966a ) . A three 
power law function was l ater introduced by [Sca lo (1986) 
[Kroupa etal.(1993)| and [Binney & Merrifield (1998). ; in 
this case the range of existence of the segments as well as 
the exponent a are functions of the investigated environment. 
After 15 years, the number of segments was raised to four 
once a universal IMF was introduced, see | Kroupa (200 1)| 
Recent developments translate the concept of IMF from our 
surroundings ( for example the first 10 pc) to open clus- 
ters such as the Pleiades, see [Moraux et al. (2004)] A con- 
nection between three segment zones and three physical pro- 
cesses has been investigated in jElmegreen (2004)] these three 
zones correspond to brown dwarf masses Ki 0.02A^q, 
to intermediate mass stars and high mass stars. The ques- 
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tion of uniformity or variability of the initial mass func- 
tion has been carefully explored, see 'Meusinger (1987)] 
IKroupa (2002)|Shadm ehri (2004)] andplmegreen (2004) 

This paper briefly reviews the physical nature of gamma- 
variate and lognormal distribution (Section l2. Il l, introduces 
an algorithm that deduces the mass from the two basic pho- 
tometric parameters B-V and My (Section ]3}, and analyses 
the IMF in the light of physical and non physical functions 
(Section]4]i. 

2. Preliminaries 

Once the histogram of the star's masses is obtained, we can 
fit it with three different probability density functions ( in the 
following pdf): two of them are well known , the lognormal 
and the power law , while the third is a gamma-variate which 
has been applied here for the first time. The lognormal dis- 
tribution /ln (a;; /J., cr) is characterised by the average value 
and the variance , a^, see [Evans et al. (2000)| The power 
law distribution has the form p{x; a) cx a;~" where a > 0. A 
refinement of the power law function consists in the multiple- 
part power law 

ax)^x-^^ , (1) 

with i varying from 3 to 4. This operation means that the 
field of existence of the function is broken into i independents 
zones. 



3 MASS DETERMINATION FROM COLOUR ABSOLUTE-MAGNITUDE 



The gamma-variate is analysed in Section l2yT1 
2.1. Fragmentation 

The stellar mass has been viewed as the Jeans mass 
in the cloud core, see for example |Lareoi^992) 



ghnegreen (1997) Elmegre en (1999)| [Chabrier (2003) 



and |Bate & Bonnell (2005)| Here conversely we adopt 
the point of view of the fragmentation that was started 
by |Kiang (1966b)| where a rigorous demonstration of the size 
distribution of random Voronoi segments was given. The 
Kiang pdf ( a gamma-variate) has the form 



H{x;c) 



r(c) 



{cxY 



(2) 



where < x < oo, c > 0, and r(c) is the gamma func- 
tion with argument c , see [Kiang (1966b)| This pdf is charac- 
terised by average value, fi=l, and variance ,a^=l/c. 

Starting from a rigorous result in ID, "Kiang (1966b )'con- 
jectured that the Voronoi cells in 2D/3D have an area/volume 
distribution represented by equation (|2j with c=4/6. 

We now briefly review a useful result that can be ob- 
tained once xi ■ ■ ■ xn , independent random variables, are in- 
troduced. Supposing that 



a;i X a;2 X 



X XN-l X XN 



N 



(3) 



it can be shown that for large N the pdf of y is approximately 
lognormal, see |Hwei Hsu (1996)| From this point of view the 
area with y = xi x X2 ot the volumes y — xi x X2 x X3 
can be fitted with a lognormal once the number 2 and 3 are 
considered a transition to large N. 

3. Mass determination from colour 
absolute-magnitude 

The (B-V) colour can be expressed as 

Tbv 



{B-V)= K 



BV 



T 



(4) 



here T is the temperature, i^BV and Tbv are two parameters 
that can be derived by implementing the least square method 
on a series of calibrated data. The bolometric correction, BC, 
can be expressed as 

Tbc 



BC = M 



bol 



Mm 



T 



10 logio T + Xbc , (5) 



where Afboi is the absolute bolometric magnitude. My is the 
absolute visual magnitude. Tec ™d A'bc two parameters 
that can be derived through the general linear least square 
method applied to a series of calibrated data. We now justify 
the theoretical basis of equations @ and We can start 
from the definition of (B-V) 



{B ~V) =niB.-mY = K - 2.5 log^ 



(6) 



where S\ is the sensitivity function in the region specified by 
the index A , K is a constant and I\ is the energy flux reaching 



the earth. We now define a sensitivity function for a pseudo- 
monochromatic BV system 



Sx^5{\-\,) i^B,V 



(7) 



where 6 denotes the Dirac delta function. On inserting the 
energy flux as given by the Planck distribution, the following 
is found 



(8) 



where Ki is a constant , c is the velocity of light , k is the 
Boltzmann's constant , h is the Planck's constant , Ab the 
wavelength of the filter B, Ay the wavelength of the filter V, 
and T the temperature. On inserting Ab = 4450 A and Ay 
= 5500 A in equation (|8j, the equation is theoretically 
derived. 

A theoretical couple, Tbv=6701 K and A'bv=-0.51 can 
be found by inserting the constants of the Planck distribution 
in equation (|8} and calibrating the relationship on the Sun 
couple, Tq =5777 K and {B — V)q)=Q.65 . In a similar way 
the bolometric correction is found to be 



BC = Mboi - A/v = if " - 2.5 log 



10 



/ SyI\dX 



(9) 



where K" is another constant. 

When the numerical value of the exponential in the 
Planck function is much greater than one, equation (|5} 
is found with the theoretical value Tbc=28402K and 
iirBC=42.45 , once BC of the sun is requested, BCq =-0.08 



atT, 







5777 K. 



The application of equation and (|5} to the calibrated 
physical parameters for stars of the various luminosity classes 
has been reported in Tabled 



Table 1. Table of coefficients from calibrated data. 

MAIN, V GIANTS, III, SUPERGIANTS I 





-0.641 ± 0.01 


-0.792 ± 0.06 


-0.749 ± 0.01 


TbvIK] 


7360 ± 66 


8527± 257 


8261 ± 67 


Kbc 


42.74 ± 0.01 


44.11 ±0.06 


42.87 ± 0.01 


rBc[K] 


31556 ±66 


36856 ± 257 


31573 ± 67 




0.062 ± 0.04 


0.32 ±0.14 


1.29 ± 0.32 




3.43 ± 0.06 


2.79 ± 0.23 


2.43 ± 0.26 



These numerical results can be visualised in Figure 
when (B-V) against T is considered and in Figure |2l 
where the relationship between BC and T is reported; 
see |Press et al. (1992)| about the fitting techniques. 

Other authors analyse the form T=T((B-V), [Fe/H], 
logj^Q g) where [Fe/H] represents the metalUcity and g 
the surface gravity, see Sekiguch i & Fukugita (2000)| or a 
two piecewise function for B-V =f (T), see for exam- 
ple Reed (1998) 

The luminosity of the star is represented through the usual 
formula 



logio(7^) = '^•4(4.74 

Lq 



(10) 



4 MASSES 




1/T [KeTvin] 



Fig.l. Least square fit of (B-V) against 1/T, MAIN SE- 
QUENCE V. Error on (B -V)=0.025. The cafibrated d ata, ex- 
tracted from Table 15.7 in' Drilling & Landolt (2000 7 are rep- 
resented with their error, and the full line reports the sug- 
gested fit. 
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Fig. 2. Relationship of EC against T, MAIN SEQUENCE V. 
Error on BC=0.025. The stars represent the calibrated data 
extracted from Table 15.7 in IDrilhng & Landolt (2000)| and 
the dotted line, the suggested fit. 



(11) 



see, for example, equation (5.120) in [Lang (1999)| The con 
nection with the mass is made through the numerical relation 
ship 

, L , , , , M 

logioll— ) = flLM + OlM logio(-r^ 

Lq Mq 

for M>Q.2Mq . 

The numerical coefficients olm and 6lm, as deduced from 
the calibrated data in Table 3.1 of |Deeming & Bowers (1984)| 
are reported in Tabled 

With the theory here adopted the mass of the star is 

logio 



M 







-0.4Mv-0.4/!:bc + 4.0 ln( 



(B-y)-ifBv 



In (10)"' 



5LM 



4 ^Bc(-(B-v)+i^Bv) ^_89g _ 

i BV 

6lm 



(12) 



From a practical point of view the percentage of reliabil- 
ity of our evaluation can also be introduced: 



e=(l- 



-) X 100, 



(13) 



where M.'^^'^ is the calibrated value of the star mass, see Table 
15.8 in Dri mng & Landolt (2000)1 and is the numeri- 

cal value here evaluated. The resulting masses are reported in 
Table|2l 

Table 2. Table of evaluated masses through an analytical 
method, equation (I12> . 



sp 


MK Class 






<%) 


BO 


V 


17.5 


17.34 


99.54 


AO 


V 


2.9 


3.42 


91.7 


FO 


V 


1.6 


1.76 


95.23 


GO 


V 


1.05 


1.13 


96.02 


KO 


V 


0.79 


0.81 


98.73 


M2 


V 


0.4 


0.38 


98.38 


K5 


III 


1.2 


6.38 


31.63 


MO 


III 


1.2 


7.15 


28.71 


AO 


I 


16.0 


22.35 


83.44 


MO 


I 


13.0 


27.73 


63.82 


M2 


I 


19.0 


28.6 


79.82 



The formulae here derived can be used to deduce the mass 
of the stars once the parameters My, B-V and luminosity 
class are provided. 

4. Masses 

We now apply the developed techniques to the nearest stars 
and to an open cluster, NGC6649. 

4.1. The masses in the first 10 pc 

The completeness in masses of the sample can be evaluated 
in the following way. The limiting apparent magnitude 
is known, for example in the case of Hipparcos it is 8. The 
corresponding absolute limiting magnitude is computed and 
inserted in equation (I12> . The limiting mass , M.^ , for stars 
belonging to the MAIN SEQUENCE V is 

logio ■ 



0.1164 m^. 



0.505 In 



0.2528 ln(0.1dj; 
0.36810" 



4.122- 



, , , 0.499 (B - V) 
,0.5107(5 - +0.3206107/ ^ ' 

where di is the maximum distance that characterises the sam- 
ple in pc. When di=\Q , m^=8 , and (B-V)=1.15 we obtain 
AA^ = 0.53A^q; this is the limit over which the sample in 
masses is complete. 

Due to the fact that Hipparcos JESA 1997b gives B-V, 
my and parallax, we can implement our algorithm on that 
database. The classification in MK classes (I,III,V) is then 
obtained by a comparison with the interval of existence of 
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4.2 The masses in open clusters 
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Fig. 3. My against (B-V) (H-R diagram) in the first 10 pc, 
this sample is incomplete in masses. 

calibrated stars. In Figure |3l we report the H-R diagram of 
the first 10 pc when the incomplete (in masses) sample is 
considered. 

In Table|3]we report the coefficients of the three pdf, aver- 
age value, the number, Ng, of stars and the merit function x^. 
A typical fit of the complete sample in the light of the power 
law pdf is reported in Figure @] conversely Figure [S] reports 
the fit through the lognormal distribution. 

Table 3. Coefficients of mass distribution in the first 10 pc , of 
the incomplete sample, GIANTS included , M > 0.2Mq, 
M = OMMq and7V,=137. 



Kiang 


x'- 


= 14.28 


c= 


1.67 


Power law 


x' 


= 17.6 




:1.85 


lognormal 


x' 


=0.95 




=0.6 o-=0.57 
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Fig. 4. Log-Log histogram of mass distribution with a super- 
position of the fitting line (the power law distribution) when 
the complete sample is considered. Parameters as in Table|3l 
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Fig. 5. Histogram (step-diagram) of mass distribution with 
a superposition of the lognormal distribution when the com- 
plete sample is considered. Parameters as in Table|3l 



Table 4. Coefficients of mass distribution in the first 10 pc, 
of the complete sample, M > 0.53Xq, M = O.98A^0, 
GIANTS excluded , and Ns=57. 



Kiang 


'2 

X = 


--12.62 


c= 


3.42 


Power law 


x' 


=8.7 


a: 


=2.42 


lognormal 


x' 


=4.55 




= 1.01 cr=0.41 



4.2. The masses in open clusters 

The distribution of masses in open clusters can be 
considered a further application of the pdf here stud- 
ied. We briefly review some of the typical clusters 
subject to investigation: IC348 ( IPreibisch et al. 20031 
; IMuench et al. 20031 ; ILuhman et al. 20031 ), 
Taurus ( ILuhman 2000l ; 'Briceno et al. 2002' ), 
Orion trapezium ( [Hillenbrand & Carpenter 200 ; 
ILuhman et al. 20001 ; IMuench et al 20021 1. PTeiades 
( IBouvier et al. 1998" ; Luhman et al. 2000> . and M35 
( p^arrado y Navascues et al. 2001, ). 



Here attention is focused on NGC6649, 
see [Walker & Laney (1987)| In this case the distance 
modulus is 11 (1585 pc) and this fact allows the determina- 
tion of the masses from Afv-(B-V) diagram, see Figure |6j 
once the algorithm outlined in Sect.jSjis adopted. 

In Table |5] we report the parameters relative to the three 
pdf. Once the framework of the three power law is accepted , 
Table |6lreports the results that are visualised in Figure0 



Table 5. Coefficients of mass distribution in NGC 6649 
where A^s=234 and JA. = 1.O4A40. 



Kiang 


x' 


=35.44 


c= 


=5.02 


Power law 


x'- 


248.49 


a-- 


=3.12 


lognormal 


x' 


=29.39 




=1.03 cr=0.32 



5 CONCLUSIONS 



Table 6. Three power law in NGC 6649 where A^,j=234. 



Mv 



1 — ' — ' — ' — ' — r- 



I 



1.1 ^^ 1.3 1.4 1.S i.« 

B-V 

Fig. 6. Mv against (B-V) (H-R diagram) in NGC6649. 




Log {JI/jHq) 



Fig. 7. Log-Log histogram of mass distribution in NGC6649 
with a superposition of the three-part power law. Parameters 
as in Table |6l 




Log {JL/JIqJi 



Fig. 8. Log-Log histogram of 100000 masses randomly gen- 
erated according to lognormal distribution (parameters as in 
Table |6} with a superposition of the three-part power law. 



Interval 




a 


0.52 < M/Mq 


<0.88 


-6.83 


0.88 < M/Mq 


< 1.8 


2.55 


1.80 < M/Mq 


<3.65 


1.0 




Log {M/ULq) 



Fig. 9. Log-Log histogram of 100000 masses randomly gen- 
erated according to the Kiang distribution with c=4 and Ml = 
2Mq. On the right , M > O.51A^0, is visible the fitting 
power law distribution with a=2.56. The vertical dashed line 
represents the limiting mass = 0.53A1q. 

5. Conclusions 

A new method represented by equation il2\ . allows us to de- 
duce the mass of stars from photometric data. This method is 
limited in mass , see equation J14l i , in the same way as the 
Hipparcos data are complete up to m„=8. 

In the second part of this paper the obtained mass distri- 
bution is fitted with the pdf of physical nature: the lognormal 
and the gamma variate. The gamma-variate is applied for the 
first time as a fitting function of mass distribution. When the 
complete sample of masses is fitted with the gamma-variate , 
c=3.42 , is found ; a comparison should be made with c=4 , 
which is the value suggested by |Kiang~1966b) in the case of 
2D fragmentation. Careful analysis of the results of Table [S] 
and Table 13 indicates that the lognormal gives the best fit to 
the star's masses. On the contrary the three power law func- 
tion simply represents the various segments of a continuous 
curve. This point of view can be expressed by plotting 100000 
stars generated according to lognormal data, see Figure|8l 

In Figure IH , the three indices ai are now varying in a 
continuous way. The great variations in visible in Table|6| 
and Figure^l can be due to the transition from a great number 
of stars ( smooth behaviour of aj ) to a small number of stars 
( jumps in a, or variability in the IMF ). 

It is now possible in the light of the limiting mass Ml^, 
see equation ( I14t , to explain the reason for which a power 
law fits the masses well in the first 10 pc. The limiting mass 
introduces a cut in the physical sample at a value that is prob- 
ably near the averaged value, and therefore only the masses 
of value greater than the averaged value are detected. This 
point of view is simulated in Figure |9] in which a real con- 
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